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This  paper  presents  an  algorithm  for  computing  the  optimal 
target  path  for  two  aircraft  traversing  a target  area  from 
different  directions.  There  are  constraints  on  the  maneuver- 
ability of  each  aircraft  which  prohibit  it  from  attacking  every 
target.  The  algorithm  chooses  a subset  of  targets  whose  de- 
struction will  yield  maximum  value  to  the  attacking  force. 

The  basis  of  the  algorithm  is  the  branch  and  bound  method, 
with  upper  bounds  computed  by  dynamic  programming.  Several 
variations  are  considered,  such  as  payload  limit,  an  increased 
number  of  aircraft  from  each  direction,  and  a three-directional 
attack.  An  example  problem  is  solved  using  the  basic  model. 

A Fortran  IV  computer  program  is  included.  Computation 
time  versus  problem  characteristics  is  discussed. 
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I . INTRODUCTION 


Optimal  assignment  of  targets  to  aircraft  on  a strike 
mission  can  greatly  increase  the  effectiveness  of  an  air  attack 
against  ground  installations.  The  purpose  of  this  paper  is  to 
present  an  algorithm  which  will,  given  the  location  and  rela- 
tive military/industrial  worth  of  key  enemy  positions,  select 
a sequenced  subset  of  targets  whose  destruction  will  yield 
maximum  value  to  the  attacking  force.  It  assumes  that  the 
decision  maker  has  complete  knowledge  of  the  targets  and  their 
value.  A stochastic  extension  can  easily  be  incorporated 
when  considering  hardened  targets.  This  is  done  by  multiply- 
ing the  probability  of  killing  a particular  target  by  its  value. 

The  initial  course  of  an  aircraft  approaching  a target 
area  is  denoted  by  $.  Due  to  anti-air  defenses,  there  is  a 
limit  placed  on  the  maximum  number  of  degrees  which  an  aircraft 
may  deviate  from  its  initial  course.  This  angular  deviation  is 
denoted  by  9.  Its  effect  is  to  restrict  the  aircraft's  movement 
at  any  point  to  a cone  whose  vertex  angle,  29,  is  bisected  by 
the  aircraft's  initial  course,  <j>.  Thus  at  any  point  (x,y)  in 
the  target  area,  the  aircraft's  course  options  are  between 
♦-0  and  $+6.  This  creates  an  ordering  among  the  targets  which 
dictates  the  sequence  in  which  they  must  be  considered.  The 
allowable  deviation  is  pictured  in  Figure  1. 
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Not  Possible 


(x,y) 


Figure  1.  Cone  of  Allowable  Course  Deviations 

Because  of  the  course  constraint,  there  are  only  a limited 
number  of  paths  through  the  target  area  ^hich  are  feasible. 

The  task  of  the  decision  maker  is  to  choose  that  set  of  tar- 
gets which  yields  the  optimal  combined  value. 

References  [1]  and  [2]  present  a method  for  determining 
the  optimal  set  of  targets  for  one  aircraft  traversing  an  area 
containing  MN  targets.  The  problem  is  formulated  as  an  MN+l 
stage  dynamic  programming  problem  [Ref.  3] . The  targets  are 
numbered  in  decreasing  order  from  an  imaginary  line  drawn  tan- 
gent to  the  boundary  at  which  the  aircraft  enters  the  target  area. 


The  stages  are  represented  by  straight  lines  drawn  through 
the  targets,  parallel  to  the  boundary  tangent  line.  Stage  n 
corresponds  to  the  parallel  line  drawn  through  target  n. 

The  starting  point  of  the  aircraft,  outside  the  target  area, 
is  stage  MN+1.  It  is  located  so  that  any  target  is  accessible 
from  stage  MN+1.  The  stage  diagram  is  depicted  in  Figure  2. 
Should  two  targets  be  equidistant  from  the  boundary  tangent 
line,  either  one  may  be  assigned  the  next  sequential  number, 
and  the  problem  reduces  to  MN  stages. 
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Figure  2.  The  Stage  Diagram  for  One  Aircraft 

Entering  from  below  the  Target  Area. 


The  state  variable  x denotes  the  lateral  position  of  the 

n 

aircraft  at  stage  n.  Dn  is  the  decision  as  to  the  heading  of 
the  aircraft  as  it  moves  from  stage  n to  stage  n-1.  Dn  is 
restricted  to  lie  in  the  set  of  feasible  headings  Sn,  from 
4>— 8 to  <j>+6.  The  state  variable  xn  is  then  a function  of  xn+1 
and  Dn+1>  This  function  is  referred  to  as  the  stage  trans- 
formation t. 

The  return  function  for  stage  n is  denoted  by  r . Letting 
Pn  be  the  lateral  position  of  target  n,  and  Vn  be  its  value, 


r (x  ) = 
n n 


if  x = p 
n n 


otherwise 


The  problem  is  then  written  as 


MN 

Maximize 

2 rn 

n-1  n n 

Subject  to: 

xn=t(xn+l'Dn+l) 

n=l , . . . , MN 

Dn  « Sn 

n=2 , . . . , MN+1 

An  efficient  algorithm  for  solving  this  problem  is  dis- 
cussed in  Section  VI,  but  the  problem  can  also  be  easily  solved 
graphically.  To  do  so,  the  course  deviation  angle  must  be 
viewed  from  the  perspective  of  the  stage  diagram.  Figure  4 
shows  whether  target  n is  feasible  for  various  values  of  xn+]_- 
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STAGE 

n+1 


Figure  4.  Cone  of  Feasibility 

Using  Figure  2 as  an  example,  the  maximum  return  possible 
for  stages  (n-1)  to  1 is  recorded  on  each  stage  line  n, 
n = 1,  ...»  MN+1,  for  every  possible  value  of  xn.  Assuming 
the  value  of  each  target  is  one.  Figure  5 gives  the  solution 
for  one  aircraft  entering  from  below  the  target  area. 

Tracing  back  from  stage  MN+1,  the  optimal  set  of  targets 
is  8,  7,  6,  4,  3,  2,  for  a value  of  six. 

Dynamic  programming  can  extend  this  problem  to  M aircraft 
by  increasing  the  number  of  state  variables  to  M [Refs.  1,  2, 
and  3] , as  long  as  all  M aircraft  are  attacking  from  the  same 
direction.  But  when,  the  aircraft  enter  the  target  area  from 
different  directions,  a pure  dynamic  programming  solution  is 
no  longer  possible,  since  the  association  between  the  target 
and  the  stage  is  no  longer  valid.  A target  in  stage  three 
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for  one  aircraft  may  be  in  stage  ten  for  another,  and  the  re 
cursive  equations  have  no  meaning. 


Figure  5.  Graphical  Solution  to  the  One-Aircraft  Problem. 


This  paper  presents  a solution  to  the  problem  of  two  air- 
craft attacking  a target  area  from  different  directions.  A 
model  is  formulated,  and  the  algorithm  developed.  An  example 
problem  is  included.  Also  included  is  a computer  program  for 
implementing  the  algorithm  and  a discussion  of  its  effective- 
ness. Because  the  terminology  becomes  quite  involved,  a 
glossary  of  terms  is  provided  in  Appendix  A. 
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II.  PROBLEM  FORMULATION 


In  formulating  this  problem,  it  is  assumed  that  if  an 
aircraft  attacks  a particular  target,  that  target  is  destroyed. 

In  other  words,  there  are  no  misses.  As  was  stated  in  Section 
I,  this  does  not  preclude  a stochastic  approach  for  hardened 
targets  where  each  has  a probability  of  being  destroyed. 

For  ease  of  notation,  the  two  aircraft  are  denoted  air- 
craft A and  aircraft  B. 

The  targets  are  numbered  as  in  Section  I,  but  now  each 
target  has  two  numbers,  one  with  respect  to  aircraft  A,  and  one 
with  respect  to  aircraft  B.  They  are  denoted  A^  and  B^  , respec- 
tively, i=  1,  ...,  MN,  j = 1,  ...,  MN.  It  is  critical  to  real- 
ize that  if  any  target  A^  has  the  same  coordinates  as  any  tar- 
get Bj , then  that  A^  is  the  identical  target  B^ . In  fact,  for 
every  A^  there  will  be  a B^  identical  to  it.  Figures  2 and  3, 
which  are  identical  target  areas,  should  clarify  this  point. 

The  value  of  A^  will  be  denoted  by  V[A^] , and  the  value  of 
B j by  V[B.]  . 

A feasible  set  of  targets  for  a single  aircraft  is  one  in 
which  the  course  required  to  go  to  each  successive  target  is 
within  the  allowable  course  deviation  cone  described  in  Sec- 
tion I.  PA (I)  denotes  one  of  the  feasible  sets  of  targets 
for  aircraft  A.  It  is  convenient  to  think  of  PA (I)  as  a 
path  of  targets  which  aircraft  A will  attack.  Any  PA (I)  com- 
pletely disregards  aircraft  B in  that  it  is  computed  as  if  only 
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one  aircraft  were  attacking.  Similarly,  PB(J)  denotes  one 
of  the  feasible  sets  of  targets  for  aircraft  B. 

The  value  of  a set  of  targets,  or  a path,  is  the  summa- 
tion of  the  individual  target  values  comprising  that  set.  The 
values  of  PA (I)  and  PB(J)  are  denoted  by  V[PA(I)]  and 
' [PB { J) ] , respectively. 

In  order  for  the  solution  of  a two-aircraft  problem  to  be 
feasible,  the  sets  of  targets  for  aircraft  A and  aircraft  B 
must  be  mutually  exclusive,  that  is,  they  must  have  no  targets 
in  common.  This  is  so  because  one  aircraft  is  sufficient  to 
destroy  the  target.  No  additional  return  would  be  realized 
by  the  other  aircraft  attacking  the  same  target.  If  this 
feasibility  constraint  were  not  required  for  optimality,  aircraft 
A (or  aircraft  B)  might  forego  the  opportunity  to  attack  other 
targets  in  order  to  attack  the  target  that  both  aircraft  have 
in  common. 

The  optimal  solution  is  achieved  when  PA (I)  and  PB(J)  are 
chosen  so  as  to 

Maximize  V[PA(I)]  + V[PB(J)] 

Subject  to:  PA(I)  r\  PB(J)  = J0. 
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III.  THE  ALGORITHM 


1 

I 


The  branch  and  bound  method  forms  the  basis  of  the  solution 
algorithm  for  the  two-aircraft  problem.  It  is  discussed  in 
theory  in  Refs.  4,  5,  6,  and  7.  Only  a small  fraction  of  the 
possible  solutions  to  the  problem  is  actually  enumerated.  The 
remaining  solutions  are  eliminated  from  consideration  through 
the  application  of  bounds  that  establish  that  such  solutions 
cannot  be  optimal. 

The  algorithm  begins  by  considering  all  possible  combina- 
tions of  paths  for  both  aircraft.  It  then  breaks  this  set  of 
all  possible  combinations  into  smaller  and  smaller  subsets,  and 
calculates  for  each  an  upper  bound  on  the  value  of  the  best 
paths  contained  therein.  The  bounds  determine  the  partitioning 
of  the  subsets  and  eventually  identify  an  optimal  path  for  both 
aircraft.  The  branch  and  bound  method  represents  the  subsets 
as  nodes  of  a tree  and  the  partitioning  of  the  subsets  as  a 
branching  of  the  tree. 

Node  one  consists  of  all  possible  combinations  of  paths  for 
both  aircraft.  Using  the  single  aircraft  dynamic  programming 
(D.P.)  method  of  Section  I,  the  optimal  path  for  aircraft  A 
is  computed  for  the  direction  from  which  A enters  the  target 
area.  This  path  is  denoted  by  PA(1) . By  the  same  method,  the 
optimal  path  for  aircraft  B is  computed  for  the  direction  from 
which  B enters.  This  resultant  path  is  denoted  by  PB(1).  It 
is  important  to  realize  that  both  PA(1)  and  PB(1)  are  computed 
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as  single  aircraft  optimizations,  and  that  all  MN  targets  are 
possible  elements  of  PB(1),  even  those  in  PA(1). 

If  PA (1)  and  PB(1)  have  no  targets  in  common,  then  the 
paths  form  the  optimal  solution  to  the  two-aircraft  problem,  the 
value  of  which  is  V[PA(1)]  + V[PB(1)],  and  the  algorithm  stops. 
However,  if  there  are  targets  in  common,  then  as  was  mentioned 
in  Section  II,  the  solution  is  not  feasible.  The  summation 
V [PA (1) ] + V[PB (1) ] represents  instead  an  upper  bound,  denoted 
by  UB(1) , for  the  optimal  solution.  The  summation  of  the 
values  of  the  points  of  intersection  of  PA(1)  and  PB(1)  is  de- 
noted by  VINT ( 1) . The  target  in  the  intersection  which  has 
the  highest  value  is  denoted  by  X^d).  This  target  could  be 
written  in  terms  of  aircraft  A or  aircraft  B.  When  the  dis** 
tinction  is  necessary,  X^fl)'  will  be  written  as  either  X^d^ 
or  XAB(1)B,  for  aircraft  A and  aircraft  B,  respectively.  Geo- 
graphically, however,  X^d),  XAB(1)A'  and  XAB(1)Bare  identical- 

The  goal  is  to  find  optimal  paths  for  aircraft  A and  B 
which  have  no  targets  in  common.  The  nature  of  the  D.P.  solution 
for  a single  aircraft  path  is  such  that  it  seeks  out  those  fea- 
sible targets  with  the  highest  value.  In  node  one,  both  A and 
B sought  XAB(1) . If  the  set  of  path  combinations  in  node  one 
was  restricted  so  that  aircraft  A had  to  take  XAB(1)A  and  air- 
craft B could  not  take  X^fljg,  the  same  target,  then  that 
point  of  intersection  would  be  eliminated.  But  perhaps  the 
optimal  solution  requires  B,  rot  A,  to  take  XA0(1).  It  could 
even  require  that  neither  path  include  XAB(1).  Therefore, 
to  include  all  possibilities,  the  set  of  node  one  is  broken 
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into  two  subsets,  one  which  requires  that  the  path  of  aircraft 


A exclude  X^d)  , giving  aircraft  B the  option  of  taking  it 
or  not,  depending  on  the  single  aircraft  D.p.  solution  for  B. 
The  other  subset  requires  that  the  path  of  aircraft  B exclude 
X.nd)#  giving  aircraft  A the  option.  Thus  node  one  branches 
to  form  nodes  two  and  three. 

The  restrictions  placed  on  each  aircraft  at  node  I are  de- 
noted by  R (I) . If,  for  example,  R(I)  = A . , B . , Bv,  the  path 
of  aircraft  A would  be  required  to  exclude  target  A^ , and  the 
path  of  aircraft  B would  be  required  to  exclude  targets  B ^ 
and  B^. 

The  branching,  with  restrictions,  is  illustrated  in 
Figure  6. 


Figure  6.  The  Start  of  a Tree 

Now  consider  node  two  with  the  restriction  vector  R(2) 
containing  the  single  element  XAB(1)A.  Again  using  the 
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singla  aircraft  D.P.  method,  only  this  time  with  the  restric- 
tion that  the  path  of  aircraft  A not  include  target  X^d)^ 
the  optimal  path  for  aircraft  A is  computed.  The  restriction 
can  be  incorporated  into  the  solution  techniques  of  Section  I 
by  temporarily  assigning  a large  negative  value  to  XAB(1), 
thereby  making  it  highly  unattractive  as  an  element  of  the 
optimal  path.  Denote  the  resultant  path  PA (2) , and  return 
the  original  value  to  XAB(1) . As  was  noted,  restricting 
PA (2)  to  exclude  XM(1)A  places  no  restrictions  on  aircraft 
B.  Therefore,  PB(2)  will  be  identical  to  PB(1).  Just  as 
with  node  one,  UB(2),  VINT(2),  and  XAfi(2)  can  be  computed. 

Node  three  is  considered  next,  with  restriction  vector 


R(3)  containing  the  single  element  X^dJg.  Since  this  means 
that  the  path  of  aircraft  B must  exclude  X^lDg,  a large 


negative  number  is  assigned  to  X._(l)_, 

Aid  B 


PB  (3) 


is  computed , 


and  the  original  value  is  returned  to  XAB(l)g.  Since  R(3) 


places  no  restrictions  on  PA(3),  it  is  identical  to  PA(1). 


Next,  UB  (3 ) , VINT(3) , and  X._(3)  are  computed. 

Ad 

A terminal  node  is  one  from  which  branching  may  still 


occur.  Nodes  two  and  three  are  terminal  nodes.  Since  only 
two  branches  may  emanate  from  any  node,  node  one  is  no  longer 
terminal,  and  need  not  be  considered  further. 


The  next  step  in  the  algorithm  is  to  choose  the  terminal 
node  which  has  the  highest  upper  bound.  Assume  it  is  node 
J (at  this  point,  J is  either  two  or  three,  whichever  has  the 


higher  upper  bound) . Node  J then  branches  to  form  nodes  four 
and  five.  The  restriction  vector  R(4)  will  equal  R(J)  with 
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The  restriction  vector 


the  addition  of  element  X._ (J) . . 

AB  A 

R(5)  will  equal  R(J)  with  the  addition  of  element  X.D(J)_. 

Ajd  d 

Again  the  restricted  paths  are  computed  and  the  algorithm 
continues  until  the  stopping  condition  is  met. 

The  branching  is  always  done  in  pairs,  a left  branch  and 
a right  branch,  as  was  shown  in  Figure  6.  No  valid  upper 
bound  comparisons  can  be  made  until  both  branchings  have 
been  performed  and  the  upper  bounds  of  both  new  nodes  have 
been  computed. 

The  stopping  condition  for  the  algorithm  occurs  when, 
following  a double  branching,  a terminal  node  is  found  whose 
upper  bound  is  greater  than  or  equal  to  all  other  terminal 
node  upper  bounds,  and  whose  paths  for  aircraft  A and  aircraft 
B are  mutually  exclusive,  that  is,  they  have  no  targets  in 
common.  This  solution  is  optimal  because  its  value  is  equal 
to  the  upper  bound  of  that  node,  making  it  the  best  solution 
possible  for  that  node.  And  since  its  value  is  at  least  as 
good  as  the  best  solution  possible  for  all  other  terminal 
nodes,  it  is  a global  optimum. 

The  steps  of  the  algorithm  are  summarized  in  Figure  7 
at  the  end  of  this  section. 

There  is  one  other  calculation  which  can  be  made  at  each 
node  which  is  of  some  interest.  The  lower  bound  for  each 
node  I is  denoted  by  LB  (I)  and  is  equal  to  UB(I)  minus  VINT  (I)  . 
This  is  a lower  bound  because  if  all  the  points  of  intersection 
on  PA (I)  and  PB(I)  were  given  to  either  aircraft  A or  air- 
craft B,  the  value  of  the  resultant  solution  would  be  LB (I). 
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If  any  terminal  node  I has  a lower  bound  on  the  optimal  value 
which  is  greater  than  or  equal  to  the  upper  bound  of  any  other 
terminal  node  J,  then  node  J may  be  completely  dropped  from 
consideration.  It  is  said  to  be  fathomed,  and  is  no  longer 
terminal.  This  will  not  speed  up  the  algorithm  or  reduce  the 
number  of  branchings  required,  since  the  algorithm  would 
never  branch  from  node  J anyway.  However,  if  computer  storage 
space  were  critical,  it  would  be  advantageous  to  incorporate 
lower  bounding,  since  once  a node  was  fathomed  it  could  be 
removed  from  storage.  It  will  not  be  employed  in  this  al- 
gorithm. 

The  algorithm  guarantees  that  an  optimal  solution  will  be 
found.  However,  it  suffers  from  a limitation  common  to  all 
branch  and  bound  methods.  For  any  untried  problem,  it  is  im- 
possible to  tell  beforehand  exactly  how  much  computation  will 
actually  be  necessary  to  find  the  optimal  solution.  Depending 
on  the  way  the  problem  is  set  up,  it  could  converge  to  the 
optimum  very  quickly,  or  for  a large,  difficult  problem  it 
could  require  such  excessive  branching  that  it  becomes  compu- 
tationally prohibitive.  This  would  be  the  case  if  the  allow- 
able course  deviation  were  very  large,  the  angle  between  at- 
tackers very  small,  and  there  were  multiple  optimal  paths. 
However,  in  sample  problems  of  target  optimization,  the  al- 
gorithm converges  very  quickly,  as  will  be  shown  in  Section 
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IV.  A TEST  PROBLEM 


Two  aircraft  are  to  attack  a target  area.  Aircraft  A 
has  an  initial  heading  of  due  north,  and  aircraft  B is  heading 
due  west  as  they  approach  the  area.  Each  aircraft  has  an 
allowable  course  deviation  of  forty-five  degrees. 

The  target  positions  and  values  are  given  in  Table  I. 

The  positions  are  given  in  terms  of  the  cartesian  plane,  with 
the  positive  Y axis  pointing  due  north. 

TABLE  I 

TARGET  POSITIONS  AND  THEIR  VALUES 


Target  Positions  (X,  Y) 


, 12 


13,  5 


, 0 


8,  6 


11,  14 


10,  11 


, 10 


Target  Values 


3,  8 

3 

1,  13 

in 

15,  15 

2 

23 
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The  stage  diagrams  for  aircraft  A and  aircraft  B,  including 
allowable  course  deviation  cones,  are  illustrated  in  Figure  8 
and  Figure  9,  respectively. 

Beginning  with  node  one,R(l)  equals  the  null  set,  meaning 
there  are  no  restrictions  on  the  path  of  either  aircraft. 

Using  the  single  aircraft  D.P.  method  for  aircraft  A,  it  is 
found  that 

PA  ( 1 ) A13  , A^^ , Ay,  A^ 

PA ( 1 ) = A^3»  A^,  Ay,  Ag  is  also  optimal,  and  either  may  be 
chosen.  For  this  example,  the  former  is  used.  Similarly, 
for  aircraft  B, 

PB  ( 1 ) — B13 , Bg , B^,  B-^ 

From  PA ( 1 ) and  PB(1),  the  following  values  are  computed: 

XAB(1)  " A3  = B1 
VINT ( 1)  = 5 

UB  ( 1)  = V [PA  (1)  ] + V [PB  (1)  ] = 29 
Since  both  paths  have  target  A3  = B-^  in  comm.cn,  they  do  not 
represent  a feasible  solution  for  the  two-aircraft  problem. 

Node  one  branches  to  form  nodes  two  and  three. 


R ( 2 ) 

- XAB(1> 

R (3 ) 

’ XAB(1> 

The  tree  at  this  point  is  illustrated  in  Figure  10. 
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Figure  10.  Start  of  the  Tree 

Considering  node  two,  PA (2)  and  PB(2)  are  computed  sub- 
ject to  the  restriction  that  PA (2)  may  not  include  target  A 
There  are  no  restrictions  on  PB(2). 

PA  ( 2 ) = Ai3'  Au»  A7  ' A6 
PB (2)  = B13,  B11#  Bg/  B5,  Bx 

The  target  of  intersection  is  Ac  = B_. 

0 O 

xab(2)  - a6  - b5 

VINT  (2)  = 5 
UB  (2)  = 29 

Moving  to  node  three,  PA (3)  and  PB(3)  are  computed  sub- 
ject to  R ( 3 ) which  states  that  PB(3)  must  exclude  target  B-, 

PA  ( 3 ) — A^3,  a^]_>  A7  i 
PB  ( 3 ) = ^11*  ® 9 ' 

2‘ 


The  target  of  intersection  is  A^  = B^. 

XAB(3)=A7=B3 

VINT (3)  = 3 

UB(3)  = 27 

Since  a double  branching  has  been  completed,  the  upper 
bounds  of  all  terminal  nodes  must  be  compared.  The  terminal 
nodes  are  nodes  two  and  three.  Node  two  has  the  highest  upper 
bound.  Since  PA (2)  and  PB(2)  have  a target  in  common,  the 
paths  cannot  be  feasible  for  the  two-aircraft  problem.  There- 
fore, node  two  branches  to  form  nodes  four  and  five.  The 
branching  tree  expands  to  Figure  11. 


Figure  11.  Expanded  Branching  Tree 

Node  four  is  subject  to  the  restrictions  from  node  two 
in  addition  to  the  restriction  that  PA (4)  cannot  include  Ag. 
PA  ( 4 ) = A13,  An,  Ag,  A5,  A4,  A2 
PB(4)  ■*  B^g , B^^,  Bg » B^,  B1 


Bg,  and 


Here  there  are  two  targets  of  intersection,  A^  = 
A2  = Bll* 

XAB(4)  " A2  " B11 

VINT  (4)  = V[A4  = Bg]  + V[A2  = Bi;l]  = 4 
UB(4)  = 28 


Node  five  has  R(5)  equal  to  R(2)  with  the  added  restric- 
tion that  PB ( 5)  not  include  xAg(2)B*  Thus  R(5)  = A3 , Bg. 

PA (5)  must  exclude  target  A^  , and  PB(5)  must  exclude  target 

B5  * 

PA(5)  » A13,  An#  A?,  Ag 
PB  ( 5 ) = B^2  , Bg  » 

There  are  no  targets  of  intersection,  but  this  does  not  mean 
that  the  optimal  solution  has  been  found,  since  terminal  node 
upper  bounds  have  not  yet  been  compared. 

xab<5)  = * 

VINT (5)  = 0 
UB(5)  = 24 


Having  completed  a double  branching,  upper  bounds  are 
now  compared  for  terminal  nodes  three,  four,  and  five,  and  it 
is  found  that  node  four  has  the  highest  upper  bound.  PA  (4) 
and  PB(4)  have  a target  in  common,  and  are  therefore  not 
feasible  for  the  two-aircraft  problem.  Since  no  other  ter- 
minal node  has  an  upper  bound  as  high  as  node  four,  the 
algorithm  branches  from  node  four  to  form  nodes  six  and 
seven. 

The  algorithm  continues  in  this  manner  with  the  upper 
bound  computed  for  each  new  node,  and  following  each  double 
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branching,  a comparison  of  terminal  node  upper  bounds  and 
a check  for  feasibility  is  performed.  When  nodes  eight  and 
nine  branch  from  node  three  and  the  bound  of  each  is  com- 
puted, it  is  found  that  node  eight  has  the  highest  upper 
bound  of  terminal  nodes  five,  six,  seven,  eight,  and  nine. 

It  is  further  found  that  PA (8)  and  PB(8)  have  no  targets 
in  common.  Therefore,  the  algorithm  stops  and  PA (8)  and 
PB(8)  form  the  optimal  solution  with  a value  of  26. 

Table  II  summarizes  the  progress  of  the  algorithm  in  the 
example.  Figure  12  illustrates  the  complete  branching  tree 
for  the  problem,  and  Figure  13  shows  the  optimal  path  of 
each  aircraft  through  the  target  area. 
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Figure  12.  Branching  Tree  for  Test  Problem 


V.  VARIATIONS 


The  model  which  has  been  formulated  can  be  modified  to 
solve  more  difficult  problems  having  additional  constraints. 

One  constraint  would  be  to  limit  the  number  of  bombs  on 
each  aircraft,  thus  limiting  the  number  of  targets  allowed 
in  the  optimal  paths  of  the  aircraft.  To  incorporate  this 
restriction,  the  single  aircraft  D.P.  method  of  Section  I 
must  be  modified.  This  can  be  done  by  increasing  the  num- 
ber of  state  variables  from  one  to  two.  At  each  stage  there 
will  be  one  state  variable,  xn,  representing  the  lateral 
position  of  the  aircraft  at  stage  n.  Another  state  variable, 
NBOMBSn,  denotes  the  number  of  bombs  remaining  in  the  air- 
craft at  stage  n.  Although  the  computation  required  for  the 
D.P.  portion  of  the  algorithm  increases  exponentially  with 
the  number  of  state  variables,  practical  problems  can  still 
be  quickly  solved  with  this  added  constraint.  The  bomb 
limitation  has  been  incorporated  into  the  computer  program 
contained  in  Appendix  B. 

In  a similar  manner,  a restriction  on  the  total  lateral 
deviation  of  the  aircraft  or  on  the  total  number  of  course 
changes  allowed  in  the  target  area  could  be  considered. 

Another  modification  would  be  to  have  M aircraft  attack- 
ing from  each  direction.  The  single  aircraft  D.P.  portion 
of  the  algorithm  becomes  instead  an  optimization  for  a single 
group  of  aircraft,  where  a separate  mutually  exclusive  path 
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is  computed  for  each  aircraft  in  the  group  (the  group  con- 
sisting of  all  aircraft  coming  from  one  direction) . The  op- 
timization for  a single  group  of  M aircraft  is  a dynamic 
programming  problem  with  M state  variables.  The  lateral 
position  of  the  ith  aircraft  at  stage  n is  denoted  by  x^n. 

The  one  directional  problem  is  discussed  in  detail  in  Ref.  2. 

To  solve  the  two  directional  problem,  the  group  of  aircraft 
coming  from  one  direction  is  viewed  as  group  A,  the  other  as 
group  B.  Since  there  will  be  no  intersection  of  targets 
within  a group,  the  only  concern  will  be  with  targets  in  com- 
mon between  the  two  groups.  As  in  Section  III,  the  highest 
valued  common  target  at  node  I,  XAfi(I),  can  be  found  and  the 
branching  performed  with  one  node  of  the  branch  restricting 
group  A to  exclude  target  X^B^A'  the  other  node  restricting 
group  B to  exclude  XAB^B‘ 

Again,  it  is  critical  to  realize  that  doubling  the  num- 
ber of  state  variables  far  more  than  doubles  the  computations 
required,  and  eventually  the  problem  will  become  computationally 
infeasible. 

The  next  modification  to  be  considered  is  the  problem 
of  aircraft  attacking  from  more  than  two  directions.  The 
general  theory  of  the  two-directicnal  problem  can  be  extended 
to  the  N directional  case,  but  the  rules  governing  branching 
become  more  involved,  and  the  number  of  nodes  required  for 
solution  greatly  increases.  One  possible  approach  where  N = 3 
will  be  briefly  considered,  with  the  aircraft  designated  A, 

B,  and  C. 
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At  each  node,  there  are  two  types  of  intersection  possible, 
a two-aircraft  intersection,  and  a three-aircraft  intersection. 
Branching  will  be  done  on  the  highest  valued  target  of  inter- 
section, whether  it  is  common  to  two  aircraft  or  three.  At 
, node  I this  point  will  be  designated  INT(I) . Assume  that  at 

node  L,  INT(L)  = A^  = B j . A double  branch  would  emanate 
from  node  L,  one  restricting  aircraft  A to  exclude  A^ , the  other 
restricting  aircraft  B to  exclude  Bj  . Should  INT(L)  instead 

equal  A.  = B.  = C,  , a triple  branching  would  be  required, 
i 3 K 

One  branch  would  restrict  aircraft  A and  B to  exclude  A.  and 

| 

Bj , respectively.  The  second  branch  would  restrict  aircraft 
A and  C to  exclude  A^  and  C^,  respectively.  The  last  branch 

I 

would  restrict  aircraft  B and  C to  exclude  B.  and  C,  , respec- 

3 * 

tively . 

The  paths  at  each  node  are  calculated  as  described  in 
Section  I,  using  the  single  aircraft  D.P.  optimization,  sub- 
ject to  the  restrictions  above.  The  upper  bound  for  any  node 
is  the  summation  of  the  values  of  the  three  paths.  The  stop- 
ping condition  is  reached,  as  in  Section  III,  when  the  terminal 
node  with  the  highest  upper  bound  has  no  targets  in  common  on 
the  paths  of  the  three  aircraft. 

If  computer  time  is  critical,  a suboptimal  solution,  as 
close  to  optimal  as  the  decision  maker  desires,  could  be  found. 
This  is  done  by  selecting  a value  which  represents  the  maximum 
difference  the  decision  maker  can  tolerate  between  the  highest 
upper  bound  and  its  corresponding  lower  bound.  When  this  value 
is  achieved,  the  suboptimal  solution  is  obtained  by  randomly 
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assigning  the  targets  of  intersection  at  that  node  to  either 
aircraft,  thus  making  the  solution  feasible. 


i 
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VI.  THE  COMPUTER  PROGRAM 

Appendix  B contains  a computer  program,  written  in  For- 
tran IV,  which  will  solve  the  target  optimization  problem  for 
two  aircraft  traversing,  from  different  directions,  a target 
area  of  up  to  100  targets.  The  difference,  a,  in  the  initial 
courses  of  the  two  aircraft,  may  vary  between  0 and  360  degrees. 
The  aircraft  may  have  between  two  and  twenty  bombs  on  board. 

The  computer  program  gives  the  user  the  number  of  branchings 
required  for  solution,  the  optimal  value  of  the  targets 
chosen,  the  path  of  targets  each  aircraft  is  to  attack,  and 
a plot  of  the  target  area  and  the  optimal  paths  through  it. 
Appendix  C contains  the  output  from  a one  hundred  target 
area,  with  an  allowable  course  deviation  of  forty-five  de- 
grees, six  bombs  per  aircraft,  and  a = ninety  degrees. 

The  input  parameters  are  the  total  number  of  targets, 
the  number  of  bombs  on  board  each  aircraft,  the  allowable 
course  deviation  angle,  the  difference  in  the  initial  courses 
of  the  aircraft,  and  the  location  and  value  of  each  target. 

The  target  positions  are  given  in  cartesian  coordinates,  with 
the  X axis  perpendicular  to  the  initial  course  of  aircraft  A. 

The  angle  ct  is  measured  counter-clockwise  from  aircraft  A, 
and  15  input  in  degrees.  The  allowable  course  deviation 
angle  is  identical  for  both  aircraft,  although  the  program 
could  easily  be  modified  to  allow  each  a separate  deviation. 

This  angle  is  also  input  in  degrees. 
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The  program  begins  by  sorting  and  numbering  the  targets, 
first  with  respect  to  the  initial  course  of  aircraft  A,  and 
then  with  respect  to  aircraft  B. 

An  MN  by  MN  matrix  is  formed  for  determining  for  either 
aircraft  whether  one  target  position  may  be  feasibly  reached 
from  another.  A "1"  indicates  feasibility,  and  a "0"  infeasi- 
bility. The  diagonal  elements  from  upper  right  to  lower  left 
are  all  zero,  indicating  that  the  aircraft  may  not  remain  at 
one  target  for  more  than  one  stage.  Denoting  any  element 
as  FEAS (I , J) , the  elements  below  the  diagonal  give  the  feasi- 
bility of  aircraft  A going  from  target  I to  target  J.  The 
elements  above  the  diagonal  give  the  feasibility  of  aircraft 
B going  from  target  J to  target  I.  This  matrix  eliminates 
the  need  to  geometrically  compute  feasibility  at  every  stage 
of  the  algorithm. 

A vector  AB  is  formed  to  correlate  the  target  numbers  with 
respect  to  A with  the  target  numbers  with  respect  to  B . If 
AB ( i)  - j,  then  A^  = B^ . 

The  path  restrictions  for  every  node  are  stored  in  a 100 
by  50  matrix  R.  The  matrix  permits  up  to  one  hundred  branch- 
ings of  the  algorithm,  and  restrictions  of  up  to  fifty  targets 
at  each  node.  Either  of  these  may  be  increased  by  the  user. 
The  restrictions  for  node  I are  stored  in  row  I of  R by  denot- 
ing a target  A^  as  negative  i,  and  a target  B,  as  j , thus 
signifying  whether  a particular  numerical  element  restricts 
aircraft  A or  aircraft  B. 
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The  single  aircraft  D.P.  portion  of  the  algorithm  is  per- 
formed in  subroutine  D.p.  of  the  program.  A sj..g’,cly  modified 
version  of  the  double  DO  loop  method  suggested  in  Ref.  2 is 
used.  It  is  presented  in  Fortran  in  a simplified  form  in 
Figure  14. 

DO  30  I = 2 , MN 
DO  20  J = 1,1-1 

C Is  it  feasible  to  go  from  target  I to  target  J? 

DO  10  K = 2 , NBOMBS 

C If  I have  K bombs  on  board  at  target  I,  is  J 
the  best  target  to  go  to? 

10  Continue 
20  Continue 
30  Continue 

C Trace  back  to  find  the  best  path 

Figure  14.  The  Simplified  Triple  DO  Loop  of 
Subroutine  D.P. 

A one-bomb  limitation  is  not  allowed,  since  the  solution 
to  the  one-bomb  problem  is  merely  to  choose  the  two  highest 
valued  targets  and  assign  one  of  them  to  aircraft  A and  one 
to  aircraft  B.  If  an  unlimited  number  of  targets  is  possible 
on  a path,  as  might  be  the  case  in  planning  a photo  recon- 
naissance mission,  the  bomb  limitation  should  be  completely 
removed  from  the  program,  rather  than  using  a very  large 
number  for  the  limitation.  This  is  so  because  increasing  the 
number  of  bombs  on  board  increases  the  computational  effort 
rorni  ired . 
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The  score  ordering  method  suggested  in  Ref.  8 was  tested 
in  the  program.  This  method  requires  that  for  target  I,  the 
targets  1-1,  ...,  1 are  stored  in  a list  in  order  of  non- 

increasing cumulated  value.  At  target  I,  the  list  is  scanned, 
starting  from  the  top,  until  a feasible  target  J is  found. 

The  cumulative  value  of  J is  then  added  to  the  value  of  I, 
and  target  I is  placed  in  the  list,  its  position  depending  on 
its  now  cumulative  value.  This  eliminates  the  need  to  scan 
all  the  lower  numbered  targets  from  target  I to  find  the 
best  one.  It  did  in  fact  result  in  reduced  computation  for  up 
to  three  bombs  on  board.  But  with  more  than  three,  the  com- 
putation required  to  update  the  list  at  each  stage  outweighed 
the  savings,  and  therefore  score  ordering  was  not  included  in 
the  program. 

For  each  node  I in  the  problem,  the  values  of  R(I),  UB(I), 
Xab ( I ) / and  VINT (I)  are  saved.  PA (I)  and  PB(I)  are  discarded 
as  soon  as  the  above  four  values  are  computed.  Since  the 
branching  tree  grows  norizontally  as  well  as  vertically  (see 
Figures  12  and  20) , a large  amount  of  storage  would  be  used 
up  in  saving  the  paths.  For  this  reason,  new  paths  are  com- 
puted at  each  new  node  from  R(I)  and  XAB(I)  of  its  predecessor 
node  I . 

Figures  15,  16,  and  17  are  time  comparisons  for  various 
input  values  of  NBOMBS , course  deviation  angles,  and  a.  The 
data  points  on  the  graphs  represent  averages  for  three  dif- 
ferent random  target  areas.  However,  the  trends  were  almost 
identical  for  each  set  of  targets. 
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Increasing  the  number  of  bombs  per  aircraft  causes  an  al- 
most linear  increase  in  execution  time.  An  increase  in  the 
course  deviation  angle  approximates  an  exponential  increase 
in  the  time  required.  Decreasing  a causes  execution  time  to 
increase.  This  should  be  expected,  since  the  closer  to  parallel 

the  two  aircraft  are,  the  larger  the  number  of  targets  in  com- 

j 

men  at  each  node  is  likely  to  be,  resulting  in  more  branching 

I 

being  required. 

4 

Figure  18  plots  execution  time  versus  the  number  of  tar- 
gets in  an  area.  Figure  19  illustrates  the  corresponding  num- 

} 

ber  of  branchings  required.  As  the  total  number  of  targets 
increases,  the  execution  times  tend  to  cycle.  One  possible 
explanation  for  this  is  that  as  the  number  of  targets  increases, 
the  amount  of  computation  in  the  dynamic  programming  subroutine 
increases.  However,  when  a path  intersection  occurs,  more 

I t 

targets  provide  more  alternative  paths  which  may  be  feasible 
and  optimal,  thereby  reducing  the  number  of  branchings  re- 
quired. For  certain  numbers  of  targets,  the  time  savings  from 
the  reduced  branchings  override  the  increased  subroutine 
computation  required,  thus  reducing  total  solution  time. 

Similar  trends  were  found  with  other  target  arrangements. 

This  again  points  out  the  unpredictability  of  the  computation 
required  for  solution. 
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NUMBER  OF  BOMBS  ON  BOARD 


TOTAL  NUMBER  OF  TARGETS  = 100 
COURSE  DEVIATION  ANGLE  = 45° 
a = 90° 


Figure  15.  Execution  Time  Versus  Bombs  Per  Aircraft 
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TOTAL  NUMBER  OF  TARGETS  = 100 
NUMBER  OF  BOMBS  PER  AIRCRAFT  = 6 
COURSE  DEVIATION  ANGLE  = 45° 

Figure  17.  Execution  Time  Versus  a 
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Figure  18.  Execution  Time  Versus  Number  of  Targets 
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Figure  19.  Required  Branchings  Versus  Number  of  Targets 


VII.  CONCLUSIONS 


The  algorithm  and  computer  program  presented  can  solve  the 
two-aircraft  target  optimization  problem  for  up  to  one  hundred 
targets,  with  an  aircraft  payload  of  up  to  twenty  bombs,  using 
a minimal  amount  of  computer  time.  In  addition,  it  can  be 
expanded  by  the  user  to  suit  his  specific  needs,  including 
more  targets,  larger  payloads,  more  aircraft,  and  a multiple 
direction  attack  scenario. 

The  algorithm  itself  could  be  improved  if  there  were  some 
way  to  recognize  the  optimal  solution  before  the  stopping 
condition  was  met.  Oftentimes,  an  early  node  will  produce  the 
optimal  path,  but  the  algorithm  continues,  because  the  node's 
upper  bound  is  not  the  highest.  This  is  the  case  where  multi- 
ple optimal  solutions  exist,  and  each  level  of  the  branching 
tree  produces  equivalent  upper  bounds.  This  is  illustrated 
in  Figure  20.  The  upper  bounds  are  indicated  above  each  node. 

At  node  13,  an  upper  bound  of  thirty-seven  is  achieved  with 
no  targets  in  common  for  PA (13)  and  PB(13).  Yet  the  branch- 
ing must  continue,  since  other  terminal  nodes  have  higher 
upper  bounds.  Finally,  by  node  25,  it  is  realized  that  the 
solution  found  at  node  13  was  in  fact  optimal.  Had  this  been 
realized  at  node  13,  the  computation  required  could  have  been 
cut  in  half. 

The  computer  program  was  developed  to  test  the  algorithm, 
and  should  not  be  considered  as  an  end-product  software  package. 


w 
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Figure  20.  A Branching  Tree  with  Alternative  Optimal 

Solutions . 

It  is  storage  inefficient,  especially  in  the  area  of  the  R ma- 
trix and  the  fact  that  information  is  kept  in  storage  for  all 
nodes  rather  than  just  the  terminal  ones.  Improvements  in  these 
areas  could  be  implemented  if  desired. 
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APPENDIX  A 


GLOSSARY  OF  TERMS 

Allowable  Course  Deviation  = the  allowable  number  of  degrees 

that  an  aircraft  is  permitted  to 

deviate  from  its  initial  course 

upon  entering  the  target  area. 

MN  = the  total  number  of  targets  in  the  target  area. 

^ Vi 

A^  = the  i target,  ordered  non-increasingly  by  its  parallel 
distance  from  the  boundary  at  which  aircraft  A enters 
the  target  area,  i = 1,  ...,  MN . 

Bj  = the  jth  target,  ordered  non-increasingly  by  its  parallel 
distance  from  the  boundary  at  which  aircraft  B enters 
the  target  area,  j = 1,  ...,  MN. 

V [A^]  = the  value  of  target  A^ . 

V [B j ] = the  value  of  target  Bj . 

PA (I)  = the  set  or  path  of  targets  for  aircraft  A computed  at 
node  I.  It  is  feasible  and  optimal  for  the  single 
aircraft  problem. 

PB  (I ) =*  the  set  or  path  of  targets  for  aircraft  B computed 

at  node  I.  It  is  feasible  and  optimal  for  the  single 
aircraft  problem. 

V [ PA ( I) ] = the  summation  of  the  values  of  the  targets  in  PA(I). 

V [PB  ( I) ] = the  summation  of  the  values  of  the  targets  in  PB(I). 

Xab(I)  = the  target  common  to  both  PA (I)  and  PB(I)  which 

has  the  highest  value. 
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XAB^A  = XAB^^  i-n  terms  aircraft  A. 

XAB ( ^ ) q = XAB  ^ ^ in  terms  °f  aircraft  B. 

VINT (I)  = the  summation  of  the  values  of  all  targets  common 
to  both  PA ( I)  and  PB(I). 

UB ( I ) = V [PA ( I ) ] + V[PB(I)]  = an  upper  bound  on  the  optimal 

solution  at  node  I. 

LB  ( I)  = UB ( I ) - VINT ( I ) = a lower  bound  on  the  optimal 
solution  at  node  I. 

R ( I)  = a set  of  restrictions  on  PA(I)  and  PB(I)  consisting 
of  a listing  of  targets  which  each  must  exclude. 

NBOMBS  = the  number  of  bombs  with  which  an  aircraft  enters 
the  target  area. 

a = the  difference  in  the  initial  courses  of  the  two  aircraft, 
measured  in  degrees  counter-clockwise  from  aircraft  A. 
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APPENDIX  B 


THE  COMPUTER  CODE 


j 

- 

j 

I I 

A TWO  DIRECTIONAL  TARGET  OPTIMIZATION  MCDEL  USING  the 
BRANCH  ANO  BOUND  ALGORITHM  WITH  DYNAMIC  PROGRAMMING. 


THIS  PROGRAM  COMPUTES  OPTIMAL  TARGET  PATHS  FOR  TwO 
AIRCRAFT  TRAVERSING  A TARGET  AREA  FROM  DIFFERENT 
DIRECTIONS.  IT  WILL  PERFORM  UP  TO  100  BRANCHINGS.  THE 
TARGET  POSITIONS  ARE  IN  TERMS  OF  THE  CARTESIAN  PLANE,  WITH 
THE  X AXIS  PERPENDICULAR  TO  ONE  OF  THE  AIRCRAFT'S  INITIAL 
COURSES.  THE  NUMBER  OF  BOM  BS  ON  BOARD  EACH  AIRCRAFT  MUST 
BE  BE’WEEN  TWO  AND  TWENTY.  THE  TARGET  AREA  MAY  HAVE  UP  TO 
ONE  HUNCRED  TARGETS. 

THE  INPUT  PARAMETERS  ARE  : 

MN  = TOTAL  NUMBER  OF  TARGETS 

NBOVBS  = NUMBER  OF  BOMBS  ON  BOAR  J EACH  AIRCRAFT 
ANGLE  = ALLOWABLE  COURSE  DEVIATION  ( IN  DEGREES  ) 

ALPHA  = DIFFERENCE  IN  INITIAL  COURSES  OF  THE  TWO  AIRCRAFT 
AX  ( I ) = X COORDINATE  OF  TARGET  I 
A Y < I ) = Y COORDINATE  Oc  TARGET  I 
V A ( I ) = VALUE  OF  T AR  GET  I 


MAIN  PROGRAM 


DIMENSION  THE  ARRAYS  AND  ASSIGN  VARIABLE  "rYP  E 

REAL  AXUOOIt  AYtlOO  ),VA(lJO)fXA<  lOCJtYAf  100)  f BX(IOO)  , 
13YQ00)  tVB(lOO)  ,X8(100)  .YBIIOO  ) ,TVA(l00  >,TVE<100), 
1VINTI100) » U8 ( 100) ,V( 100) 

INTEGER  FEAS  (100, 100  ),  R (100t  50)  , AS  ( 100)  , XAB(  ICC) , 

1 OATH (1 00)  ,P  A( 100  ) 

INTEGER  COUNT 


READ  IN  the  NUMBER  OF  TARGETS,  ThE  NUMBER  OF  3QMBS  PER 
AIRCRAFT,  THE  DEVIATION  ANGLE,  AND  TnE  CIFFERENCE  IN 
INITIAL  COURSES  0*  THE  TWO  AIRCRAFT. 

READ  (5,10)  mm,NBOMBS,  ANGLE, ALPHA 
10  FORMAT  ( 2 14 ,2  F5 .0  ) 


READ  IN  THE  TARGET  POSITIONS  AND  VAlJES 
CO  2 0 I =1  , MN 

RE  AC  (5  ,30)  AX  ( I ) , AY  ( I ) ,V  A l I ) 

20  CJNTINUE 
30  FORMAT  (3F12.0) 
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SORT  TARGETS  FOR  AIRCRAFT  A 
MNM  1 = MN-1 
DO  5 0 I=1,MNM1 
IP1  = 1+1 
IQPT  = 0 
CO  AO  J=IP1,MN 
IF  (AY(  I I .GE.AY(  Jll  GO  TO  40 
TEMPY  = A Y(  I ) 

TE M FX  = AX  (I  ) 

TEMPV  = V A 1 1 ) 

AY ( I)  = AY  ( J ) 

AX<  I)  = AX(  J) 

VA  ( I ) = VA(J) 

AY  ( J ) = T EMPY 
AX  ( J)  = TEMPX 
VA ( J > = TEMP  V 
40  CONTINUE 
50  CONTINUE 


SORT  TARGETS  C0R  AIRCRAFT  3 

ALPf-A  = ALPHA*3.  1415*326  54/180.0 
DO  60  I=1,MN 

3X  ( I ) = AX(  I ) *COS<  ALPHA)+AY(I  )*SI  N(  ALPHA) 
8Y  (I)  = -AX  ( I )*S  INI  ALPHA)+AY(  I )*C0S(  ALPHA  ) 
V8  ( I ) = V A(  I ) 

60  CONTINUE 

CO  30  I=1,MNM1 

IP  1 = I+i 

DO  70  J=IP1,MN 

IF  ( BY  ( I ) .GE  . BY  ( J ) ) GO  TO  70 

TEMPY  = 3 Y ( I ) 

TEMPX  = 3X11) 

TEMFV  = V B ( I ) 

BY(I)  = 9 Y(  J ) 

3X  (I  ) = 3X  (J  ) 

VB(  I)  = V 3 ( J > 

3Y  ( J ) = TEMOY 
ex  ( J ) = T EMPX 
V B ( J ) = TEM  P V 
70  CONTINUE 
80  CONTINUE 


COMPUTE  FEAS  I 31  LI TY  MATRIX 

ANGLE  = < (90.  C-ANGLE)*3. 141  592654/130. 0)-1.0E-3 
FEAS  (1,1  ) = 0 
CO  120  1=2, MN 
I M 1 = 1-1 
FEAS ( I, I ) = 0 

oo  no  o=i, imi 

IF  ( AX ( I ) . E9  . AX( J I ) GO  TO  90 
IF  ( BX  ( I ) . EQ . 8X  ( J ) ) GO  TO  100 

S4.0PEA  = ATAN(ABS(  ( AY  ( J)-AY  (I  ) ) / ( AX  ( J)-AX  ( I ) ) ) ) 
SLOPES  = AT  AN  ( ABS  ( ( BY  ( J )-B  Y (I))/(BX(J)-3X(  I)))) 
FE  AS ( I , J)  = 0 
CE  AS  ( J , I ) = 0 

IF  (SLOPEA.GT  .ANGLE)  FEAS(I,J)  = 1 
IF  (SLOPEB.GT  .ANGLE)  F=AS(J,I)=i 
GO  TO  no 
90  FEAS ( I , J ) = 1 
GO  TO  110 
100  FE A S ( J , I ) = 1 
110  CONTINUE 
120  CONTINUE 
C 
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C SET  UP  AB  CORRELATION  VECTOR 
DO  140  L=1,MN 
DO  130  M=1,MN 
TWOPIE  = 2.0*3.141552654 
TP  MA  = TW  CPIE— ALPHA 

1=  <ABS(AX(L)-(8X(M)*C0S(TPMA)+3Y(M)*SIN(TPMA))).GT. 
11 . 0E-2 ) GO  TO  130 

IF  ( ABS ( AY  < L ) - (-BX  < M )*  S I N (T  °M  A ) +BY  ( M )*C0S  I T PM  A ))  > .GT  . 
11.0E-2)  GO  TO  130 
ABIL)  = M 
GO  TO  140 
130  CONTINUE 
140  CONTINUE 


BEGIN  BRANCH  AND  BOJND 


ZERO  OUT  RESTRICTION  MATRIX 
DO  160  L= 1, IOC 
CO  150  M=l,50 

R ( L * M ) = 0 

150  CONTINUE 
160  C0NTI NUE 


DETERMINE  R VECTOR  FOR  NODE  I 
XAB ( 1 ) = 0 
I = 0 
J = 1 
K = 0 

170  I = 1*1 

IF  ( I . EO . 101)  GO  TO  490 
DO  1 30  L=  1, 50 
LAST  = L 

IF  ( R I 0 » L ) .EQ .0  > GO  TO  150 
IF  ( L .EQ  .50  ) GO  TO  510 
R(  I fL>  = R(  J,L) 

ISO  CONTINUE 

190  IF  (I  .EG. 2 ) K = 0 

Ic  (K.EQ.l)  GO  TO  20 0 
R(I,LAST)  = - XAB  ( J ) 

C-0  ~ TO  210 

200  R<I  ,LAST)  = A E (X  A B ( J ) ) 

K = Q 


ADJUST  TARGET  VALUES  TO  EXCLUDE  TARGETS  IN  R(I> 
210  00  230  L=  1,  50 

IF  ( R ( I , L ) . E Q .0 ) GO  TO  240 
IF  { R ( I f L ) .GT.O ) GO  TO  220 
""V  A ( — R ( I f L ) ) = VA(-R<I,L)  ) 

VA  ( — R ( I , L ) ) = -1000C.0 
GO  TO  23  0 

220  TVB  ( R ( I t L ) ) = V8(R(  I,U  ) 

V3( R(  I, L ) ) = -10000. 0 
230  CONTINUE 


CALL  DP  FOR  A AND  B 
240  IC  HK  = 1 

CALL  DP  (FEASfVA,  ICHK,MN,NBOMBS,  PATH, W, COUNT) 

VP  A = W 

NP  A = COUNT 

DO  250  L= 1 , N ° A 

PAIL)  = P A"H  ( L ) 

250  CONTINUE 
IChK  = 2 

CALL  DP  ( FEASfVB,ICHK,MN  fiNBOMBS  , PATF,W, COUNT  ) 
I F ( IOPT  .EO  .1  ) GO  TO  350 
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C RESET  TARGET  VALUES 
DO  270  L=1 » 50 

IF  (R(I,L ).EQ.0>  GO  TO  280 
IF  ( R ( I » L ) . GT  .0  ) GO  TO  260 
VA ( -R ( I ,L)I  = TV  A ( — R ( I ,L>  ) 
GO  TO  270 

260  V3 ( R ( I * L)  i = TVB<  R(  I,L)  ) 
270  CONTINUE 


COMPUTE  XABII),  VINT(I),  AND  UB  ( I ) 

280  XMAX  = 0.0 
SUM  = 0.0 
DO  3 CO  L = 1»  NP  A 
CO  290  M=I, COUNT 

IF  ( AB(PA (L ) ) .NE. PATH! M) ) GO  TO  290 
SUM  = SUM  +V  A ( P A ( L ) ) 

IF  (VA(PAIL)  ) .LE.XMAX)  GO  TO  300 
XMA  X = VA  (PAIL)  ) 

X A E ( I ) = P A ( L ) 

GO  TC  300 
290  CONTINUE 
300  CONTINUE 

VI  NT(  I ) = SUM 

UB ( I ) = VPA+W 

IF  (I.EQ.  l.AND.VINT(I).EQ.O.O)  GO  TO  340 
IF  (K.rQ. 1)  GO  TO  1 7 C 
UB  ( J ) = 0.0 


CHOOSE  FIGHEST  UPPER  BOUND  OF  ALL  TERMINAL  NODES 
ANJ  CHECK  IF  IT  IS  FEASIBLE 
UBMAX  = 0 .0 
00  310  L=  It  I- 

IF  ( UB C L)  .L5. UBMAX)  GO  TO  310 
UBMAX  = UB  (L  ) 

>i  = L 

310  CONTINUE 

IF  (V  INT  { J ) . EG .0 . 0 ) GO  TO  330 


CHECK  FCR  TIE  FOR  FIGHEST  UPPER  BOUND. 
IF  THERE  IS  A TIE.  IS  IT  FEASIBLE? 

CO  320  L*l»  I 

IF  (UB(L)  . N E . U B ( J ) ) GO  TO  320 
IF  ( VINTI  L> . NE.O. 0)  GO  TO  320 
J = L 
GO  TC  330 
320  CONTINUE 
GO  TO  170 


IF  NECESSARY.  RECOMPUTE  OP T I M AL  PATH 
33  0 NODES  = I 

IF  ( I .EQ  .J  ) GO  TO  350 
I = J 
I3PT  = 1 
GO  TO  210 
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C PRINT  0(jT  THE  NUMBER  OF  BRANCHING  NODES,  THE  OPTIMAL 
C VALUE,  AND  THE  OPTIMAL  PATHS. 

340  NODES  = I 
350  WRITE  <6,360) 

36G  FORMAT  ( ' 1 ' ) 

WR  ITE  ( 6,370)  NODE  S 

370  FOR  NAT  (•  OPTIMAL  SOLUTION  FOUND  IN', 13,'  BRANCHINGS', 
1 ////////) 

WR  ITE  <6,  380  ) UB ( J ) 

380  FORMAT  ('  THE  VALUE  OF  THE  OPTIMAL  SOLUTION  IS  *,F12.4 
1, ////////) 

WRITE  (6,  390) 

390  FORMAT  ('  OPTIMAL  AIRCRAFT  PATHS',///) 

WRITE  (6,  400) 

400  FORMAT  (6X.«  AIRCRAFT  A',///) 

WRITE  (6,410) 

410  FORMAT  ( 4 X , • X • ,12X , • Y • , // ) 

DO  420  1=1, NPA 
XA  ( I ) = AX(  PA  (I  ) ) 

Y A ( I ) = AY  ( P A ( I ) ) 

WRITE  (6,430)  XA  ( I ),YA(  I) 

420  CONTINUE 

430  FORMAT  (1  X,  F3  .2,  5X  , F8 .2  , / ) 

WRI  TE  ( 6,440  ) 

440  FORMAT  (///, 6X, 'AIRCRAFT  3',///) 

WRITE  (6,410) 

DO  450  1=1, COUNT 

XB  ( I ) = BX  ( P ATH  ( I ) ) *C3S  (TPM  A)  +8  Y ( PATH(  I ) ) * SIM  ( TPM  A) 
Y3(I)  = -BX(  PATH<  I ) )*SIN(TPMA  I EY  (P  ATH  ( I ) ) *CCS(  T^WA  ) 
WRITE  ( 6,  430)  XB(  I ) , Y9( I ) 

450  CONTINUE 


PLOT  THE  TARGETS  AND  THE  PATHS 
WRITE  (6,460) 

460  FORMAT  ( » 1'  ) 

CALI  PLOT  p (AX, AY,  MN, II 
CALL  PLOT  P (XA,YA,NPA,2 > 

CALL  PLOTP  (X  e,  Y 9,  COUNT  ,3  ) 

WRITE  ( 6,470) 

470  FORMAT  ( ///,40X,  'PATH  OF  AIRCRAFT  A DENOTED  BY  •*■') 
WRITE  (6,480) 

48C  FORMAT  (/,40X,'PATH  OF  AIRCRAFT  3 DENOTED  BY  *', 

1 ////////////) 

GO  TO  530 


CHECK  FOR  DIMENSIONAL  ERRORS 
490  WR  ITE  ( 6,500  I 

500  FORMAT  (////,'  REQUIRED  BRANCHING  EXCEEDS  IDO  NODES', 
1////) 

GO  TO  530 
510  wRITE  (6,520) 

520  FORMAT  (////,«  RESTRICTION  VECTOR  EXCEEDS  50', ////I 
530  STOP 
END 


56 


onn  non  on  o o<“>o  <~>oo  oo 


1 


SUBROUTINE  OP  ( F E AS , V , ICHK , MN ,NBOMB S , PATH , W, C CUNT) 

THIS  SUBROUTINE  COMPUTES  THE  OPTI m*l  PATH  FOR  A SINGLE 
AIRCRAFT  TRAVERSING  THE  TARGET  AREA. 


REAL  VV(  100,20  >,  V(  100)  ,CVAL  ( 20) 

INTEGER  P (100 ,20) ,FEAS (100,100), MOVE (20 ) , PATH (100 ) 

INTEGER  C CUNT 

CVAL(l)  = O'.O 

MOVE ( 1 ) = 200 

DO  20  K=1,NB0MBS 

VV (1, K)  = V (1) 

P( 1 ,K)  = 200 
20  CONTINUE 


I IS  THE  TARGET  THAT  THE  AIRCRAFT  IS  AT  PRESENTLY 
00  70  I =2  ,MN 
VV<  I, 1)  = V(  I ) 

CO  30  K=2,NBOMBS 
CVAL(K)  = 0.0 
MJVE(K)  = 200 
VV  (I,K)  = V ( I) 

30  CONTINUE 
IM1  = 1-1 


J IS  THE  TARGET  THAT  THE  AIRCRAFT  MAY  GO  tq  FROM  I 
00  50  J -1  , 1 Ml 

1=  ( ICHK.5Q.1  .AND.  FEAS(  I,  J)  .EQ.O)  GC  TO  50 
IF  (ICHK.E0.2  .ANO.FEASI  J,  I ) .EQ.O  ) GO  TO  50 


K IS  THE  NUMBER  OF  BOMBS  LEFT  AT  TARGET  I 
CO  AO  K=2  » NB  OMBS 

IF  ( ( VV(  I ,K ) + VV( J , ( K-i)  ) ) .LE.CVALIK  ) ) GO  TO  AC 
CVAL(K)  = VV  ( I,K  )+VV  (J  , (K-l  ) ) 

MOVE  ( K)  = J 
AO  CONTINUE 
50  CONTINUE 

00  60  K = 1 , NB  OMBS 

IF  (CVALU  ) .NE.0.0  ) VV  ( I » K ) =C  VAL  ( K ) 

P(  I , K ) = MCVE<K) 

60  CONTINUE 
70  CONTINUE 


COMPUTE  FIRST  TARGET 
CVAL  (NBOM  BS)  =0.0 
DO  80  J=1,«N 

IF  ( VV ( J , NBCMBS ) . LE • CV  AL ( NBOM  BS  ) ) GOTO  30 
CVAL(NBOMBS)  = VV(J,NB0M8S) 

MOVE  (NBOM  BS  ) = J 
80  CONTINUE 

* = CVAL(  NB0M8S) 


COMPUTE  OPTIMAL  PATH 

PATF(l)  = MOV  E (NBOM  ES  ) 

COUNT  = 1 
LEFT  = NB0M3S+1 
CO  90  I =2  , M N 
COUNT  = I 
LEFT  = LEFT-1 

PATH(I)  = P( ( FATH(I-l)  ),L  EFT l 
Ic  (P( (PATH( I )),( LEFT-1  )) .EQ. 200)  RETURN 
90  CONTINUE 
RETURN 
ENO 
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THE  COMPUTER  OUTPUT 


OPTIMAL  SOLUTION  FQUNO  IN  7 BRANCHINGS 


THE  VALLE  OF  THE  CPTIMAl  SCILTICA  IS  216.2231 


OPTIMAL  AIRCRAFT  PATHS 

AIRCRAFT  A 


X Y 


A 7 . 7 1 

48.91 

31.71 

99.13 

171.61 

297.32 

317. C5 

445.32 

639.03 

886.22 

754.96 

985.23 

AIRCRAFT 

e 

X 

Y 

794.25 

411.59 

461. 19 

389.44 

424.49 

262.81 

363.66 

378.05 

328.99 

324.02 

209.62 

363.97 
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PATH  OF  AIPCRAFT  A DENOTED  BY 
PATH  CP  A I PC  BAFT  B DENOTEC  BY 
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